On the Determination of the Rotational Oblateness of Achernar 



A. C. Carciofi ^, A. Domiciano de Souza^, A. M. Magalhaes^, J. E. Bjorkman^, F. Vakili^ 

ABSTRACT 

The recent interferometric study of Achernar, leading to the conclusion that 
its geometrical oblateness cannot be explained in the Roche approximation, has 
stirred substantial interest in the community, in view of its potential impact in 
many fields of stellar astrophysics. It is the purpose of this paper to reinterpret 
the interferometric observations with a fast rotating, gravity darkened central 
star surrounded by a small equatorial disk, whose presence is consistent with 
contemporaneous spectroscopic data. We find that we can only fit the avail- 
able data assuming a critically rotating central star. We identified two different 
disk models that simultaneously fit the spectroscopic, polarimetric, and inter- 
ferometric observational constraints: a tenuous disk in hydrostatic equilibrium 
(i.e., with small scaleheight) and a smaller, scaleheight enhanced disk. We be- 
lieve that these relatively small disks correspond to the transition region between 
the photosphere and the circumstellar environment, and that they are probably 
perturbed by some photospheric mechanism. The study of this interface between 
photosphere and circumstellar disk for near-critical rotators is crucial to our un- 
derstanding of the Be phenomenon, and the mass and angular momentum loss of 
stars in general. This work shows that it is nowadays possible to directly study 
this transition region from simultaneous multi-technique observations. 

Subject headings: polarization — interferometry — spectroscopy — stars: emis- 
sion line. Be — stars: individual (Achernar) 
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Introduction 



Interferometry greatly i ncrea sed our knowledge of the Be Star Achernar [a Eri, HD10144). 
Domiciano de Souza et al.l (l2003l . h ereafter D03) observ ed Achernar during fall 2002 with 
the ESO-VLTI/ VINCI instrument jKervella et aliboosh in the K band and found that the 
star is highly oblate. By converting the individual visibilities into equivalent uniform disk 
angular diameters they d erived a ratio between the maxi r num a nd minor elongation axis 
a/b = 1.56 ± 0.05. Lat er. iKervella fc Domiciano de Souzal (120061) detected a tenuous polar 
wind and very recently Kervella fc Domiciano de Souzal ( 2007 ) found that Achernar has a 
main- sequence, lower mass companion. 

The available interferometric data, together with the abundance of data from the lit- 
erature, forms a body of information that is probably unrivaled by that of any other Be 
star. There is, however, much that is still unknown about the star and its circumstellar 
environment. As pointed out by DOS, the determination of the actual stellar oblateness from 
the interferometric data is not a simple task, since there are so many unknowns involved. 
In their original work DOS, based on simultaneous spectroscopic data, assumed that there 
was no circumstellar material at the time of the observations, and that the observed shape 
was purely photospheric in origin. Their modeling of the observations with a rotationally 
deformed, gravity darkened star, assuming the Roche approximation of uniform rotation and 
centrally condensed mass, led to the conclusion that the stellar flattening required to explain 
the observations exceeds that of a critical rigid rotator . 



Recent theoretical developments in the theory of differentially rotating stars (e.g.. |Jackson et al. 



20041 ) show that they can have equatorial radii larger tha n the Roche lin ii t, wh ich is physically 
compatible with the interferometric results. However, IVinicius et al.l (120061 ) demonstrated 
the existence of a tiny, yet non negligible emission in Ha at the time of the VLTI observa- 
tions, an indication that there was some circumste llar material around Achernar, probably 
associated with a small rotating disk. Furthermore, IVinicius et al.l (120061 ) showed, from tem- 
poral analysis of profile variations of the line He I A6678 A, that orbiting gas clouds are a 
very frequent feature of Achernar, even during its quiescent phase. 

The presence of this circumstellar material raises the question of how much it can 
contribute to the observed visibilities. In this paper we investigate how the presence of a 
small disk around Achernar may alter to the interferometric signal, and the effects this may 
have on determination of fundamental parameters such as stellar rotation rate and stellar 
flattening. 
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2. Models 

To compute the emerg e nt spectrum of the system we use the computer code HDUST 
dCarciofi fc Biorkman|[2006l : ICarciofi et aPbooe h. This code solves the coupled problem of 
the non-local thermodynamic equilibrium (NLTE) and radiative equilibrium for arbitrary 
three-dimensional envelope geometries, using the Monte Carlo method. Given a prescription 
for the star (rotation rate, shape, flux distribution, effective temperature, luminosity, etc.) 
and for the circumstellar gas (density and velocity distribution) the code calculates the NLTE 
Hydrogen level populations and the electron kinetic temperature. Once those quantities are 
determined, the emergent spectrum (spectral energy distribution, line profiles, polarization, 
synthetic images, etc.) is calculated. 

For the present work, we modified HDUST to treat non-spherical, gravity darkened 
stars. We assume the stellar shape to be an spheroid, whi ch is a reasonable approximation 
for the shape of a rigidly rotating star (IFremat et al.ll2005l ). Given the stellar rotation rate, 
Q/Qcrit, the ratio between the equatorial and polar radii is determine d from the Roche 
approximation for the stellar surface equipotentials (IFremat et al.ll2005r). F or the gravity 



von ZeipeJ 1924 ). according to 



darkening we use the standard von Zeipel flux distribution (j^ 
which F{9) oc ges{0) oc T^q{9), where ges and Tes are the effective gravity and temperature 
at stellar latitude 6. In HDUST, the star is divided into a number of latitude bins (typically 
50) which have an associated T^s and Qp fj and emit wi th an spectral shape given by the 
corresponding Kurucz model atmosphere (iKurucz Ill994j ). 



For the disk density distribution we assume the following expression 



eq 



■cq 



ZU 



exp {-z^/H'^) 



(1) 



which is similar to the density distribution of an is othermal viscous decretion disk in vertical 
hydrostatic equilibrium (e.g. ICarciofi et al.l 120061 ). In the above equation, w is the radial 
distance in cylindrical coordinates, i?eq is the equatorial radius of the star and po is the disk 
density scale. We write the disk vertical scaleheight, if, as 



J. 5 



■eq; 
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where Hq is the scaleheight at the base of the disk. For isothermal disks in vertical hydrostatic 
equilibrium Hq 

{GM/R, 



a/2 



aV^j.^lReq, where a is the sound speed and the critical velocity. Vent 



■cqj 



is the Keplerian orbital speed at the stellar surface. 



We chose a set of fixed stellar parameters, listed in Table [H The value for the polar 
radius came directly from the interferometry, once we fixed the stellar inclination angle to 
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be 65deg after the recent work of ICarciofi et al.l (120071 ) . It is worthy noting, however, that 
the results we show below are little affected if we allow the inclination angle to vary within a 
reasonable range (say, ±5deg). The value for the stellar luminosity was obtained by fitting 
the available photometric data (see caption of Figured]). 

Each model of the star plus disk system has 4 free parameters: (1) the stellar angular 
rotation rate, Q/Qcrit, (2) the disk density scale, po, (3) the disk outer radius, R^, and (4) 
the disk scaleheight, Hq. 

We have two well-known observational constraints, viz., the 2002 interferometric obser- 
vations (DOS) and the 2002 Ha line profile (jVinicius et al.ll2006l . fig. 12). Unfortunately, no 
contemporaneous measurement of the linear polarization exists, but we can impose, with a 
reasonable level of conf idence, an upper limit for the polarization based on the recent results 



of ICarciofi et al.l (120071 ). In that paper, the results of a polarization monitoring of Achernar, 



carried out between 2006 July and November, is reported. Those measurements were taken 
during a period in which A chernar was active, w ith a tenuous circumstellar disk, as indicated 
by the weak Ha emission ( ICarciofi et al.l 120071 . fig. 3). W e adopt as an u pper limit for the 
polarization in our modeling the lowest value reported in ICarciofi et al.l (120071 ). which was 
0.12% in the B band. As we shall see below, fixing this upper limit for the polarization level 
has important consequences for the modeling. 



3. Results 

Our modeling procedure is as follows. For a given set of Q/Qcnt, Po, Rd and Hq we 
calculate synthetic images in the K band, centered at A = 2.15 /im, the Ha line profile, the 
continuum polarization and the spectral energy distribution (SED). The model amplitude 
visibilities where obtained from the Fourier transform moduli of the model images at 2.15 fim, 
normalized by the total K band flux of the synthetic image. 

Since the Ha emission from October 2002 is very small and the true photospheric profile 
of Acher nar is not known, we model, instead of the observed profile, the residual emission 



profile of IVinicius et al.l (120061 ) . fig. 12. This profile was obtained by subtracting from the 
October 2002 observations the average profile of the 1999 period, which is believed to be 
purely photospheric. The equivalent width (hereafter EW) of this emission profile is —0.29 

A. 

In Figure [2] we show the visibility curves along the polar and equatorial directions 
together with the corresponding K band images for five representative models that were 
chosen to illustrate different aspects of our solution. In Table |5] we list the model parameters. 
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along with some model results, such as the B band polarization level, the Ha EW and the 
K band flux excess, E^, which is deflned as F^/F^ — 1, where F^ is the stellar flux without 
the disk at wavelength A. 

Let us begin discussing models 1 to 3, for which the star was assumed to be rotating 
critically. We have used n/n^nt = 0.999999 (Kq/Krit = 0.9993) instead of 1 to avoid the 
unphysical situation of having Teg = in the stellar equator. Model 1 corresponds to a small 
and relatively dense disk, with Hq given by Eq. [21 i.e., model 1 corresponds to a disk in 
vertical hydrostatic equilibrium. For this model the visibility curves matches very well the 
observations and the polarization is within the adopted limit, but the Ha emission is very 
weak, as a result of the small disk size. 

One way to increase the Ha emission is to make the disk larger. For model 2 we adopted 
a larger value for R^, but the density had to be decreased in order to keep the polarization 
within the adopted limit. This model reproduces well both the visibility curves and the Ha 
EW. 

Another way of increasing the line emission is to raise the disk scaleheight. If one raises 
if by a factor of, say, k, the disk density scale must be decreased by a factor k^^^ to keep the 
polarization approximately constant. Model 3 corresponds to our best scaleheight enhanced 
model for a critically rotating star. The visibility curves and the Ha EW are well reproduced 
by the model, and the polarization is within the adopted upper limit. 

We also studied models with lower values of Q/Qcrit- In this case, our attempts to flt 
simultaneously the visibihty curves, polarization and Ha EW were unsuccessful. 

Models 4 and 5 are examples of the results we have obtained for sub-critical stars. 
Model 4 corresponds to a small and dense disk in vertical hydrostatic equilibrium. This 
model reproduces quite well the visibilities, but have too large B band polarization and 
no Ha emission. A scaleheight enhanced model with lower density, such as model 5, have 
better values for polarization and Ha EW, but does not have enough IR flux in the equator 
to account for the observed size. 

Our lack of success on fltting the observations with sub-critical models is a direct result 
of two physical effects, both related to the presence of gravity darkening. For the critical 
models 1 to 3, the bolometric stellar flux at the equator is about 1100 times lower than the 
flux at the pole. Since the disks are relatively small and the star oblate, light from the pole 
cannot reach the disk and, as a result, the net scattered flux is small. As Q/Q^rit decreases, 
the flux at the equator rises, and the star becomes less oblate. Both mechanisms cause a 
signiflcant increase in the scattered flux and, thus, in the polarization level. 



- 6 - 



4. Discussion 

The presence of a small disk around Achernar is a realistic possibility to explain the 
interferometric observations. As shown by DOS (see also Figured]), Roche models without a 
disk cannot account for the observed aspect ratio. 

The combined constraints imposed by the interferometry, spectroscopy and the adopted 
upper limit for the polarization have allowed us to impose very narrow limits on the model 
parameters, within our assumption of a rigidly rotating star in the Roche approximation 
surrounded by a small disk. The main result is that the star must be rotating very close to 
critical, since all of our models with n/r2crit < 0.992 that successfully reproduced both the 
visibilities and the Ha EW (e.g. model 5) had too large polarization levels, for the reasons 
discussed above. 

We have found two critical models that reproduce equally well the observations: a large 
and dense disk in hydrostatic equilibrium (i.e., geometrically thin; model 2) and a small and 
more tenuous disk with enhanced scaleheight (i.e., geometrically thick; model 3). Because the 
lack of contemporaneous polarization measurements, the model parameters shown in Table [2] 
have a degree of uncertainty that is difficult to estimate. This stresses the importance of 
having simultaneous multi-technique data. 

The current paradigm for Be star disks is that of a geometrically thin viscous Keplerian 
disk in vertical hy drostatic equilibrium. This paradigm ha s been corroborated by several 



recent studies (e.g. ICarciofi et al.ll2006l : iMeilland et al.l 120071 ) 



When a star is close to critical rotation, the stellar photospheric material near the 
equator, which is weakly bound to the star because of the strong centrifugal force, can 
eventually escape the star provided it is given an extra energy by some other mechanism 
(e.g. stellar pulsation, interaction in a binary system, photospheric activity). This material 
is likely to have, initially, a very complicated density and velocity distribution. As the gas 
diffuses outwards as a result of viscosity, the density and velocity distributions will tend to 
relax and become in hydrostatic equilibrium. 

Thus, it is reasonable to assume that there is a transition region between the photosphere 
(which is characterized by large densities of the order of 10~^° g cm^) and the disk itself, 
which has densities at least 10 times lower. Many unknowns exist concerning the size and 
physical properties (e.g. density, temperature, velocity field) of this smooth transition region. 

With the available data we were able to narrow down substantially the range of possible 
values for Achernar's model parameters, but we cannot distinguish yet between model 2 and 
3, i.e., between a model in hydrostatic equilibrium and a model with enhanced scaleheights 
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that might result from the perturbation of the gas by some photospheric mechanism. 

As shown in Figure [H high-precision spectropolarimetry might add invaluable informa- 
tion, since models 2 and 3 have different values for the size of Balmer and Pashen jumps 
in polarization. Also, with simultaneous spectro-interferometry and spectropolarimetry one 
can spatially resolve the velocity field, thus establishing whether the material is in equilib- 
rium or not. The observed residual emission profile (Figure [1]) is indicative of velocities fields 
more complex than the simple Keplerian rotation that was assumed in this work. 

In any case, one important consequence of this work is that we have shown that with 
current observing techniques and state-of-the-art modeling it is already possible to study 
the properties of the very inner layers of the disk in Achernar and possibly other nearby Be 
stars. Future simultaneous spectropolarimetry, spectro-interferometry and photometry will 
allow us to study the properties of this region and determine, for instance, the size of the 
transition region and the point at which the disk becomes in vertical hydrostatic equilibrium. 

This work was supported by FAPESP grant 04/07707-3 to ACC, FAPESP grant 01/12598- 
1 to AMM and NSF grant AST-0307686 to the University of Toledo (JEB). AMM also 
acknowledges partial support by CNPq. 
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Table 1. Fixed Stellar Parameters. 



Parameter 


Value 


Reference 


-Rpole 


7.3 Rq 


This work 


Luminosity 


3150 Lq 


This work 




350 km s~^ 


Vinicius et al. f2006) 


inclination 


65 deg 


Carciofi et al. r2007) 


distance 


44.1 pc 


Perrvman et al. fl997l 



Table 2. Model Parameters. 



Model 


^/^crit 


Po (g cm- 


-3) 


Rd 


HiRe) 


Rp 


Ie. 


Ha EW (A) 


Pb (%) 


Ek 


1 


0.999999 


1.0 X 10- 


-11 


13.7 


0.45 


1.5 


5.8 


0.10 


0.02 


0.17 


2 


0.999999 


3.8 X 10- 


-12 


19.1 


0.45 


1.5 


5.8 


-0.34 


0.10 


0.14 


3 


0.999999 


2.6 X 10- 


-12 


16.1 


1.1 


1.5 


5.8 


-0.27 


0.12 


0.12 


4 


0.992 


1.0 X 10- 


-11 


13.7 


0.45 


1.4 


1.8 


0.20 


0.31 


0.15 


5 


0.992 


1.7 X 10- 


-12 


13.7 


0.91 


1.4 


1.8 


-0.05 


0.15 


0.03 
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Fig. 1. — Emergent spectrum for models 1 to 5. Top left: visi ble SEP. The red tr i- 
angles correspond to BVR photometry from the NOMAD catalog (jZacharias et al.ll2005l ). 
for which no observational errors were available. T op right: IR SED. The JHK photom- 
etry is from the 2MASS catalog (jCutri et al.l 120031 ) and the mid-IR photometry is from 
Kervella fc Domiciano de Souza ( 20071 ). The observational errors for the mid-IR data is of 



the order of 2%. Bottom left: continuum polarization. Bott om right: Ha e missio n profile. 
The red curve corresponds to the residual emission profile of IVinicius et al.l (120061 ). 
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Fig. 2. — Squared visibilities for the models of Table [2] (solid lines), along the polar (upper 
curves) and equatorial (lower curves) directions. The dashed lines represent the K band 
visibilities for uniform disk (UD) angular diameters of 1.62 mas and 2.53 mas derived by 
DOS; the corresponding ±l-sigma uncertainties are shown as light-color bands. These UD 
diameters indicate the maximum (equatorial direction) and minimum (polar direction) sizes 
measured on Achernar in the K band with VLTI/VINCI. The vertical dotted lines indicate 
the maximum baseline available from the VLTI data. The insets show the model image in 
the K band, in logarithmic scale. The bottom right plot shows the squared visibilities for a 
gravity darkened Roche star at near-critical velocity (models 1 to 3) without the disk, and 
indicates, as shown by DOS, that such model does not reproduce the observations. 



